Virtual covariance matrix reconstruction-based adaptive beamforming for small aperture array

Recently, many robust adaptive beamforming (RAB) algorithms have been proposed to improve beamforming performance when model mismatches occur. For a uniform linear array, a larger aperture array can obtain higher array gain for beamforming when the inter-sensor spacing is fixed. However, only the small aperture array can be used in the equipment limited by platform installation space, significantly weakening beamforming output performance. This paper proposes two beamforming methods for improving beamforming output in small aperture sensor arrays. The first method employs an integration algorithm that combines angular sector and gradient vector search to reconstruct the interference covariance matrix (ICM). Then, the interference-plus-noise covariance matrix (INCM) is reconstructed combined with the estimated noise power. The INCM and ICM are used to optimize the desired signal steering vector (SV) by solving a quadratically constrained quadratic programming (QCQP) problem. The second method proposes a beamforming algorithm based on a virtual extended array to increase the degree of freedom of the beamformer. First, the virtual conjugated array element is designed based on the structural characteristics of a uniform linear array, and received data at the virtual array element are obtained using a linear prediction method. Then, the extended INCM is reconstructed, and the desired signal SV is optimized using an algorithm similar to the actual array. The simulation results demonstrate the effectiveness of the proposed methods under different conditions.


Introduction
Adaptive beamforming technology has also been known as spatial filtering anti-interference technology, and there have been numerous research breakthroughs in this field in recent years [1][2][3].The principle of this technology is to process sample data received by each sensor according to certain beamforming criteria and using various algorithms.After performing the weighted superposition of data received by each array element, the main lobe of the beam is aligned with the signal of interest (SOI), and a null is formed in the direction of the interference signal to suppress the interference, which effectively improves the system's output performance.Because most methods based on this technology require prior information on the incoming direction of the desired signal, they have been frequently combined with the high-resolution direction of arrival (DOA) estimation algorithms [4,5] to improve an antenna array system's interference suppression ability and output performance.A minimum variance distortionless response (MVDR) beamformer can maximize the output signal-to-interferenceplus-noise ratio (SINR) under the ideal condition that the SOI steering vector (SV) and interference-plus-noise covariance matrix (INCM) are known [6,7].However, the effect of adaptive beamforming can significantly decrease in many practical applications compared with the ideal situation due to the presence of SOI in the sample covariance matrix (SCM), array geometry errors, and DOA mismatch.
Recently, several robust adaptive beamforming (RAB) techniques have been proposed to reduce the effect of model mismatches and improve the robustness of beamforming.For instance, the eigenspace-based method [8][9][10][11] can help to improve the robustness of beamforming algorithms.However, noise pollution at a low signal-to-noise ratio (SNR) can degrade performance to a certain extent; also, the determination result of the subspace dimension can affect the output of this method.Diagonal loading (DL) [12][13][14] has been a common method for increasing the robustness of beamforming.However, the setting of the diagonal loading level determines the performance of beamformers.Although the adaptive diagonal loading algorithm [15,16] can obtain a loading level adaptively, it is challenging to maintain robustness in a wide range of input SNR.In [17][18][19][20], the uncertainty set constraint technique was proposed to attain higher performance.The aforementioned methods are effective at improving the output SINR under certain error conditions.
Further, to minimize the desired signal's impact on adaptive beamformers, the INCM reconstruction method [21,22] was proposed to improve output performance by weakening the SOI in the SCM.In [21], the authors first proposed the idea of INCM reconstruction, where the INCM was reconstructed by using the Capon power spectrum integrated over an interference-plus-noise sector region, and the SV of the desired signal was optimized by solving a quadratically constrained quadratic programming (QCQP) problem.In [23], a low-complexity INCM reconstruction method based on spatial power spectrum sampling was developed.The method in [24] applied maximum entropy power spectrum to reconstruct the INCM by integrating over the angular sector of interference-plus-noise as well as the desired signal region.In [25], an effective RAB method was developed by reconstructing the INCM based on a combination of power method processing and spatial spectrum matching.In [26], the distributed digital subarray antennas were used to form a contiguous virtual array after the gaps between the subarrays were filled with virtual array elements.Then, the INCM of the contiguous virtual array was reconstructed.[27] proposed an extended INCM reconstruction method based on linear prediction to generate virtual sensor data and extend array aperture.In addition, combining the INCM reconstruction method and the coprime array [28,29] can further enhance the output performance of beamforming.
What's more, the array aperture size also affects the beamforming performance.Under particular experimental conditions, such as small detection or communication devices, the spatial aperture of an array can be smaller than the wavelengths involved [30].For linear arrays, the array aperture size is related to both the interval of elements and the number of elements.When the number of array elements is fixed while the interval of elements decreases, or when the interval of elements is fixed while the number of elements decreases, the array aperture will decrease, so the beamforming performance will decline.In order to improve the output performance of beamforming algorithms on small aperture arrays, this paper proposes two RAB approaches based on the INCM-reconstructed method.In the first proposed method, a RAB algorithm based on an actual array is proposed.This algorithm combines the angular sector integration method and line searching integration method along the gradient vector to reconstruct the interference covariance matrix (ICM), and the INCM is constructed by combining the obtained ICM with the estimated noise covariance matrix.Then, the reconstructed ICM and INCM are used to optimize the desired signal SV by solving the QCQP problem.In the second method, a RAB algorithm based on a virtual extended array is proposed.First, the virtual conjugated element is generated according to the structural characteristics of a uniform linear array, and data received at the virtual array element are predicted using the linear prediction method.Further, the INCM of the extended array is reconstructed using the ICM reconstruction algorithm in the first method, and the desired signal SV is optimized.The main contributions of this paper include: 1: Two beamforming algorithms for small aperture arrays are proposed.The first proposed method combines the angular sector integration method and line searching integration method along the gradient vector to reconstruct the ICM, which achieves well performance in different conditions.
2: The second proposed method utilizes the linear prediction method to get a virtual extended array.Based on the ICM reconstruction and the desired signal SV optimization algorithm in the first method, the virtual element is used to expand the array's aperture and improve the degree of freedom of the array.
The remainder of this paper is organized as follows.Section II describes the signal model and provides necessary background information about adaptive beamforming.Section III introduces the proposed RAB methods.Section IV presents simulation results.Finally, Section V concludes this paper.

Signal model
Consider a uniform linear array (ULA) with M elements that receive Q far-field uncorrelated narrowband signals, including one desired signal and Q − 1 interference signals.Then, the kth array received signal can be expressed as follows: where the waveform of signal s 1 is regarded as the SOI, and a 1 is the corresponding SV; s q and a q denote the waveform and SV of the qth interference signal, respectively; n is the additive white Gaussian noise.
An MVDR beamformer obtains an optimal weighting vector by minimizing interference and noise power without distorting the SOI, which can be expressed as follows: where R in is the INCM, which is given by: where R i and I M represent the ICM and identity matrix, respectively; s 2 q and s 2 n are the powers of the qth interference signal and white Gaussian noise, respectively.
The solution to the optimal weighting vector is given by: Because the actual INCM and SV are unavailable in real applications, the INCM can be replaced in two ways.One is to use the SCM R, which can be written as follows: where K is the number of snapshots.
Another way is to reconstruct an INCM based on the available information.
The output of the beamformer is given by: In the remainder of this paper, the hypothetical DOA of the SOI and qth interference are respectively denoted by ŷ1 and ŷq , where q = 2, 3, . .., Q.The desired signal SV obtained using array geometry and hypothetical DOA information is denoted by â1 , and the qth interference SV is denoted by âq .

Proposed methods
Generally, the gain of an array system is proportional to the array aperture size.As the array aperture size decreases, the array system's gain will also decrease.The output performance of small aperture arrays can be improved by designing more robust beamformers or extending arrays with virtual elements.

Beamforming algorithm based on actual array
The covariance matrix reconstruction method, which uses a spatial power spectrum integral over a particular sector as the estimate of the matrix, can effectively improve the robustness of beamforming.Inspired by this, the prposed1 method reconstructs the ICM over an angular sector and a gradient vector set.The estimated noise power is calculated based on the relationship between residual noise and actual noise.Based on the obtained INCM, the nominal SV of the desired signal can be optimized by solving a QCQP problem.ICM reconstruction.In addition to the method defined by Eq (5), a real INCM can also be replaced by reconstructing the INCM.The integral algorithm in the angle region of a signal has a strong robustness to the direction error and can reconstruct the ICM by integrating the Capon space spectrum in the angular sector of the interference, which can be expressed as follows: where Θ q represents the angular sector of the qth interference signal, and PðyÞ is the Capon spectrum in the θ direction, which can be expressed by [31]: Then, the ICM reconstruction in Θ q can be approximately performed using a discrete sum method as follows: where θ l , l = 1, 2, . .., L is the sample points in Θ q , and L represents the number of discrete sample points.
The basic idea of reconstructing the ICM by Eq ( 9) is to use the Capon power spectrum integral in the angular sector where the interference is located.The Capon power spectrum is mainly distributed in the actual signal SV [32].In addition, it can be considered that multiple virtual interference signals are added near the interference estimation direction ŷq .Although this reconstruction method can overcome the effect of direction estimation error on the calculation result to a certain extent, the acquisition of PðyÞ and âðyÞ requires accurate array calibration.If other array errors exist, the ICM reconstruction by Eq (9) will lead to not ideal results.The algorithm presented in [32] adopts a volume integral method for the INCM reconstruction to reduce the sensitivity to array calibration error.However, the number of discretization points in the integration region of [32] is two to the power of M. When the number of matrix elements is large, the number of discretization points will be too large to calculate.
Further, to address other errors and decrease the number of discretization points in the integration region, the ICM reconstruction along the gradient direction from âq is proposed.The basic idea of this method is that there can be multiple contour vectors with the same Capon power value near an actual signal SV, and âq is one of them.Therefore, if several nominal SVs are selected in the small range near ŷq to form a hyperplane, the estimated SV âq can be regarded as a tangent point between the hyperplane and the contour line.Then, from the tangent point âq , the true SV can be found along the gradient direction perpendicular to the hyperplane [33].Based on this inference, the ICM can be constructed from the tangent point âq in the gradient direction using the line integral method.
Next, the gradient can be obtained by using the subspace-based method to reconstruct the scan matrix of the SV âq in a certain range of directions as follows: where [θ 1 , θ 2 , . .., θ M−1 ] is a set of (M − 1) sample points uniformly distributed in Θ q .
Then the non-full rank matrix H q H H q can be constructed, and its eigendecomposition is performed as follows: where γ q1 � γ q2 � . . .� γ qM represents the eigenvalues, and b qm is the eigenvector corresponding to γ qm . Let , which is a hyperplane space containing M − 1 eigenvectors corresponding to the maximum of (M − 1) eigenvalues.Since matrix H q H H q is not full rank, the eigenvectors corresponding to the first (M − 1) large eigenvalues are sufficient to span the space where the hyperplane is located.According to the property of the conjugate symmetric matrix, the eigenvector b qM is orthogonal to the other eigenvectors.Further, b qM is orthogonal to the hyperplane space, so it can be used as a gradient vector estimation of âq [33].
The estimated nominal interference SV is adjusted by adding a scaled gradient vector: Inspired by [33], the gradient-based reconstruction ICM can be expressed as a discrete sum of the adjusted SV, which can be expressed as follows: Combining Eqs ( 9) and ( 13), the final ICM can be reconstructed as follows: Noise power estimation and SV optimization.Using the idea presented in [34], the noise power for the noise covariance matrix can be estimated as follows: where � s 2 n is the residual noise power, which is computed by: where θ j is a discrete sample point in the noise region Θ n , and J is the number of sample points.
Then, the reconstructed INCM is given by: The desired signal's SV can be optimized by adopting the idea of maximizing the beamformer's output power [21].In addition, to prevent the estimated desired SV from converging to other interference regions, the ICMs obtained in two different ways, by Eqs ( 9) and ( 13), are both considered when designing the constrain conditions: where v s ? is one of the components of v s , and it is orthogonal to â1 ; v s is the mismatch vector of â1 , and it can be decomposed into two components, namely v s jj and v s ?; v s jj is parallel to â1 and does not affect the beamforming quality.Therefore, v s jj is not considered in Eq (18).
The optimization problem in Eq (18) represents a quadratically constrained quadratic programming (QCQP) problem and can be solved by a convex optimization toolbox.The optimized SV is expressed as � a 1 ¼ â1 þ v s ?, and the weight vector obtained based on the actual array can be calculated by:

Beamforming algorithm based on virtual extended array
By using the real data received from an actual array, a virtual array can construct the received data at the virtual element to expand the array.The proposed2 method extends the ICM reconstruction method in proposed1 to the virtual extension array using the virtual element.Compared with the original actual array, the advantage of a virtual extended array is that it can increase the number of degrees of freedom.
Extended array data generation.A virtual extended array based on the MVDR method has three tasks to solve.The first task is to determine the expansion mode of a virtual array, that is, to determine the expansion array's SV.The second task is to generate the virtual array's received data, and the third task is to construct the INCM of the extended array.To extend the aperture and degree of freedom of an array, this study designs a virtual extended array using the conjugate symmetry of ULA [35].The extended array is shown in Fig 1.
If the first element x 1 of a real array is the point of symmetry, then the real and virtual elements are axially symmetric about x 1 .The number of sensor elements in an extended array is (2M − 1), of which (M − 1) elements are virtual elements.Assume that the signal SV of an actual array with DOA θ is âðyÞ ¼ ½1; e À j2pd sin y=l ; . . .; e À j2pðMÀ 1Þd sin y=l � T ; then, the SV of a virtual element is estimated to be expressed from right to left as follows: âV ðyÞ ¼ ½e j2pd sin y=l ; . . .; e j2pðMÀ 1Þd sin y=l � T ð20Þ where λ represents the signal wavelength, and d is the spacing between adjacent sensors.According to Eq (20), the SV of a virtual element is conjugate of the SV of an actual array.Thus, the extended-array SV estimation can be expressed as follows: ãðyÞ ¼ ½e j2pðMÀ 1Þd sin y=l ; . . .; e j2pd sin y=l ; 1; . . .; e À j2pðMÀ 1Þd sin y=l � T ð21Þ For the purpose of facilitating the subsequent calculation, the received data of a virtual extended array are expressed as follows: xðkÞ ¼ ½x À ðMÀ 2Þ ðkÞ; x À ðMÀ 3Þ ðkÞ; . . .; x 0 ðkÞ; x 1 ðkÞ; x 2 ðkÞ; . . .
Generally, except a DOA estimation algorithm can predict the DOA, information on each signal in an array is unknown, so the received data at a virtual element cannot be calculated directly by Eq (1).However, a linear prediction method can be used to generate a virtual array element's data [27,36].For a conjugate virtual extended array, the mth forward prediction sensors can be estimated in an iterative way as follows: where m = 0, 1, 2, . .., M − 2, α ¼ RÀ 1 F d F is the forward prediction coefficient, and RF is the cross-correlation matrix of the input matrix for a linear system, which is expressed by: In the forward prediction method, the first real element is usually regarded as reference data, and the second to Mth elements are considered the input matrix, which is expressed as Further, d F is the cross-correlation between x F (k) and x 1 (k), and it is defined by: Extended INCM reconstruction and SV optimization.The extended SCM can be expressed as follows: where R is a semi-positive definite matrix with a dimension of (2M − 1) × (2M − 1), and its condition number is defined by: where εmax and εmin are the largest and the smallest eigenvalues of R, respectively.It should be noted that when the number of virtual sensor array elements increases, the dimension and condition number of R will increase accordingly.A high condition number causes numerical instability in matrix inversion, which can lead to errors in Capon spectrum estimation [27].
Since the condition number is related to the matrix's eigenvalue, the eigenvalue can be changed using a diagonal loading algorithm to reduce the condition number.Assume that RL is the covariance matrix modified by applying a diagonal loading method on R; then, it holds that RL ¼ R þ iI 2MÀ 1 .Therefore, a new condition number resulting from the diagonal loading level ι can be expressed as follows: To avoid relying on experience when designing ι, cond v+ι can be limited by the condition number cond r of R. Further, define cond v+ι � cond r , such that: where εmax and εmin are the maximum and minimum eigenvalues of R, respectively.The value of a diagonal loading level needs to satisfy the following condition: Therefore, the extended Capon spectrum in θ can be expressed as follows: According to Eq (7), the extended ICM related to the incident direction of the qth interference can be reconstructed as follows: Similarly, the gradient-based extended ICM corresponding to the qth interference can be reconstructed as follows: where � ãq ¼ ãq þ Z bqM , and bqM is the gradient vector estimation obtained based on the qth interference signal's extended SV.
Combining Eqs ( 32) and ( 33), the final extended ICM can be reconstructed by: Furthermore, the estimated extended INCM can be calculated by: In this study, the extended SV of the desired signal is optimized using the method defined by Eq (18) as follows: where Rd ¼ P Q q¼2 Rdq , and Re ¼ P Q q¼2 Req ; ã1 is the extended SV estimation of the desired signal; ṽs? is one of the components of ṽs , and it is orthogonal to ã1 ; ṽs is the mismatch vector of ã1 .The convex optimization toolbox can also be used to solve the QCQP problem in Eq (36).The extended optimized SV is expressed as � a 1 ¼ ã1 þ ṽs ?, and the weight vector obtained based on the virtual extended array can be calculated by: Finally, the corresponding extended array beamforming output is given by:

Simulation results
In this study, a ULA with M = 10 omnidirectional sensors was used.There were four independent signals impinging from the directions of 5˚, −30˚, 30˚, and 50˚.It was supposed the first  [33], INCM reconstruction method with subspace decomposition, steering vector estimation and correction [38] (INCM-SDC), the extended virtual method based on linear prediction [27] (LPV), and signal vector estimation [26] (SFV).In this section, the proposed beamforming method that uses Eq (19) to obtain the weight vector is referred to as proposed1, and the method that uses Eq (37) to determine the weight vector is denoted as proposed2.
The interference and SOI angular sector were Y q ¼ ½ ŷq À 5 � ; ŷq þ 5 � � and , respectively; the complement sector of interference and SOI was defined as a noise region Θ n .A uniform sample angular interval 0.1˚was used in all angular sectors.The number of SVs constructed for H q in the actual experiment was the same as the number of samples in Θ q to achieve better performance.The number of sampling points in Θ int was I = 40 for the INCM-volume method.The number of dominant eigenvectors was set as seven in the INCM-subspace method.The parameter ξ in the INCM-MCP and proposed methods was 0.1 with an interval of 0.01.The number of virtual elements in the LPV was set to 20.The virtual element in the SFV was obtained by using the virtual conjugate extended array.
In the SINR calculation, the INCM of all virtual array methods was obtained using the data received by the real extended array.When the SNR of the SOI changed, the number of snapshots was fixed to K = 1, 000.When the number of snapshots varied, the SNR was fixed to 10 dB.The sample data used in all simulations contained the SOI at all times, and all the results are the average of 200 Monte Carlo runs.The QCQP problem was solved using the CVX [39].

Exactly known signal steering vector
This example considered the situation where the exact steering vector was known.Fig 2 displays the output SINR of the proposed methods versus the SNR, and Fig 3 shows the deviations between all tested methods and the optimal SINR.The performance of the SFV beamformer based on a virtual extended array was higher than the proposed algorithms and optimal values in this test.This is because there is no error in the SV, and when the information of the SV is completely known, this method can obtain more accurate incident signals.Then the ICM obtained by combining the incident signals with the known extended SV will be close to the ICM of the real extended array.The output of the extended virtual array beamformer LPV was also higher than the optimal values when the SNR was less than 6 dB.However, the output performance of LPV deteriorates as the input SNR increases.The output SINR of the tested beamformers versus the number of snapshots is presented in Fig 4, where it can be seen that the output SINR results of the SFV and LPV were better than those of the proposed methods and other tested beamformers.insensitive to the look direction error in one incident direction, while the look direction error in the other directions affected its output to a certain extent.

Mismatch caused by multiple model errors
In this case, look direction error, gain error, phase error, and location perturbation were all considered.The look direction error was assumed to be uniformly distributed in [−2˚, 2˚].The gain error and phase error were randomly distributed according to N(0, 0.05 2 ) and N(0, 0.1π 2 ), respectively, and the sensor location error was uniformly distributed in the interval [−0.05, 0.05] measured in sensor space.The mean output SINR versus the input SNR for a fixed number of snapshots is given in Fig 6 .Fig 7 shows the deviation of each beamformer from the optimal beamformer.This shows that the proposed1 method performs well in the whole SNR range while the performance of the proposed2 method decreases when the input SNR increases.The output SINR versus the number of snapshots for fixed input SNR is given in Fig 8 .The results for the proposed methods were close and better than the other techniques.

Mismatch caused by coherent local scattering
In this scenario, the effect of a mismatch in the desired signal SV caused by coherent local scattering [40] on the array output was examined.Under this condition, the actual desired signal SV was expressed as a ¼ aðy 1 Þ þ P 4 p¼1 e jφ p aðy Lp Þ, where a(θ 1 ) is the direct path, a(θ Lp ) is the pth coherently scattered path in a direction θ Lp , and all θ Lp , p = 1, 2, 3, 4 are randomly distributed in a Gaussian distribution with mean of 5˚and standard deviation of 1˚; φ p represents the path phase, whose value was independently and randomly generated from the inter-

Performance comparison for different arrays apertures
This experiment was conducted to analyze the performance of beamformers at different array apertures.The input SNR was fixed at 10 dB, the number of snapshots was set to 1, 000, and the input signals were unaltered.The look direction mismatch corresponding to the signals was uniformly distributed in the range of [−2˚, 2˚].First, the sensor number was fixed at 10, and it was examined how the test methods performed when the inter-sensor spacing was altered.The output SINR results of different methods versus the ratio of inter-sensor spacing d to signal wavelength λ are presented in Fig 12, where d/λ represents the ratio of array spacing to wavelength.The proposed1 and proposed2 methods could achieve stable performance when d/λ was higher than 0.3.While the output of the INCM-volume method and the LPV algorithm tended to be stable when d/λ was larger than 0.4, and the LPV algorithm was stable only when d/λ was from 0.5 to 0.7.The output SINR results of the other methods fluctuated with the value of d/λ.
Next, the interval of elements was set to half wavelength to verify the output performance of the test algorithms for a different number of elements.The output SINR results of all tested algorithms versus the number of elements are presented in Fig 13, where it can be seen that the output results of the INCM-volume, proposed1, and proposed2 methods increased with the number of components.The proposed algorithms had the best performance among all methods when the number of array elements was less than nine.

Performance comparison under extra degrees of freedom
The last example analyzed the beamforming performance under extra degrees of freedom (DOF).In this test, the number of input signals increased to 11, the sensor number was fixed to 10, and the spacing between array elements was half wavelength.The information on each input signal is shown in Table 1.It was supposed that the first signal was the desired signal.
Since the LPV and proposed2 method are processed on the virtual extended array, both can realize the beamforming when the number of incident signals exceeds the number of array elements.The output SINR results of different methods versus the input SNR and the number of snapshots are shown in Figs 14 and 15, respectively, where it can be seen that the output of the proposed2 method was better than that of the LPV method.The LPV output results were stable when the number of snapshots was larger than 300.In addition, the output of the proposed2 method was better when the number of snapshots was less than 400 than when it was higher than 400.

Conclusions
This paper proposes two RAB methods based on covariance matrix reconstruction denoted by proposed1 and proposed2.The proposed1 method performs angular sector integration and gradient vector integration to reconstruct the ICM and combines the estimated noise covariance matrix to reconstruct the INCM.Then the rebuilt ICM and INCM are used to optimize the desired signal's SV.Simulation results demonstrate that this method can achieve good performance under different errors and a small array aperture by improving the robustness of the  beamformers.The proposed2 method employs a virtual extended array and uses the conjugate property of a linear array.Data received at a virtual array element are produced by the linear prediction method, and the INCM of the extended array is reconstructed based on the ICM reconstruction algorithm similar to that used in the proposed1 method.The proposed2 method performs well for different array apertures and can increase the number of degrees of freedom of the beamformer by adding virtual elements.The complexity of the proposd2 method is higher than that of proposd1 due to the introduction of the extended array.To ensure beamformer robustness and consider computational efficiency, the proposed1 method should be considered when the number of signals is less than the number of array sensors.When the number of signals exceeds the array sensors, the proposed2 method should be considered.However, the output of the proposed2 method decreases as the look direction error increases.As a result, the robustness of the proposed2 method is an important issue to investigate in our future work.

Fig 8 .
Fig 8.The output SINR versus the number of snapshots.https://doi.org/10.1371/journal.pone.0293012.g008 val of [0, 2π].The look direction mismatch corresponding to the signals was uniformly distributed in the range of [−1˚, 1˚].It should be noted that the look direction mismatch, θ Lp , and φ p changed from run to run but were fixed from snapshot to snapshot.The output SINR results of all tested methods versus the input SNR are presented in Fig 9, and the deviations between them and the optimum result are presented in Fig 10.The results demonstrated that the proposed1 had the best performance among all methods, especially for the SNR values from −10 dB to 10 dB.This mismatch did not affect the results of the